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ABSTRACT 

From extensive radiative transfer calculations we find that clumpy torus models with Mo ~ 5-15 dusty clouds 
along radial equatorial rays successfully explain AGN infrared observations. The dust has standard Galactic 
composition, with individual cloud optical depth Ty ~ 30-100 at visual. The models naturally explain the 
observed behavior of the 10/im silicate feature, in particular the lack of deep absorption features in AGN 
of any type. The weak 10/im emission feature tentatively detected in type 2 QSO can be reproduced if in 
these sources jVq drops to ~2 or 7y exceeds ^100. The clouds angular distribution must have a soft-edge, 
e.g., Gaussian profile, the radial distribution should decrease as l/r or 1/r 2 . Compact tori can explain all 
observations, in agreement with the recent interferometric evidence that the ratio of the torus outer to inner 
radius is perhaps as small as ~5-10. Clumpy torus models can produce nearly isotropic IR emission together 
with highly anisotropic obscuration, as required by observations. In contrast with strict variants of unification 
schemes where the viewing-angle uniquely determines the classification of an AGN into type 1 or 2, dumpiness 
implies that it is only a probabilistic effect; a source can display type 1 properties even from directions close 
to the equatorial plane. The fraction of obscured sources depends not only on the torus angular thickness but 
also on the cloud number Mo- The observed decrease of this fraction at increasing luminosity can be explained 
with a decrease of either torus angular thickness or cloud number, but only the latter option explains also the 
possible emergence of a 10/im emission feature in QS02. X-ray obscuration, too, has a probabilistic nature. 
Resulting from both dusty and dust-free clouds, X-ray attenuation might be dominated by the dust-free clouds, 
giving rise to the observed type 1 QSO that are X-ray obscured. Observations indicate that the obscuring torus 
and the broad line region form a seamless distribution of clouds, with the transition between the two regimes 
caused by dust sublimation. Torus clouds may have been detected in the outflow component of FLO maser 
emission from two AGN. Proper motion measurements of the outflow masers, especially in Circinus, are a 
promising method for probing the morphology and kinematics of torus clouds. 

Subject headings: dust, extinction — galaxies: active — galaxies: Seyfert — infrared: general — quasars: 
general — radiative transfer 



1. INTRODUCTION 

Recent VLTI inter ferometric observation s in the 8-13 fj,m 
wavelength range by Tristram et al. (2007) confirm the pres- 
ence of a geometrically thick, torus-like dust distribution in 
the nucleus of Circinus, as required by unification schemes 
of Seyfert galaxies. Several aspects of their data require that 
this torus is i rregular, or clumpy, in agree ment with the earlier 
prediction of iKrolik & Begelmarj d!988l) . 

We have recently developed the first formalis m for handling 
clumpy AGN tori and presented initial results dNenkova et al.1 
I2002t lElitzur et al.ll2004t lElitzur||2006l 120071) . The reported 
clumpy models have since been employed in a number of 
observational studies, including the first analy sis of Spitzer 
obser vations by the GOODS Legacy project (Treis teret al.l 
2004). Our clumpy torus models were also employed 
in the analysis of spatially-resolved, near-diffra ction-limited 
10_um spectra of the NGC 1068 nucleus (Mas on et al.1 
120061). The geometry and kinematics of both water maser 
dGreenhill & Gw inn 1997; Gallimore et al. 2001) and narrow- 
line emission (ICrenshaw & Kraemerl 120001) indicate that the 
NGC 1068 torus and a ccretion disk are oriented nearly edge- 
on. The Maso n~et~ai1 clumpy model for IR emission is the 
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first to correctly reproduce the observed near-IR flux with an 
edge-on orientation. In contrast, smooth-density models re- 
quire viewing angles 22°-30° above the equatorial plane in 
order to bring into view the warm face of t he torus backsid e 
dGranatoetal.lll997t iGratadour et afj|2003t iFritz et alj|2006l) . 
dumpiness is also essential for understanding the puzzling 
interferometry result that dust temperatures as different as > 
800 K and ~ 200-300 K are found at such close proximity 
to each other dSchartma nn et alj 12005). The mounting ob- 
servational evidence in favor of clumpy, rather than smooth, 
dust distrib ution in AGN tori has sparked add itio nal modeling 
efforts bv lDullemond & v an Bemmel (2005) and lHonig et alj 
(120061) . 

This two-p a per series expands the analysis of 
i Nenkova et alJ (12002b . In its first part dNenkova et alJ 
I2008L part I hereafter) we develop the full formalism for 
continuum emission from clumpy media and construct the 
source functions of dusty clouds — the building blocks of the 
AGN torus. Here we assemble these clouds into complete 
models of the torus, and study the model predictions and 
their implications to IR observations. In comparing the 
predictions of any torus model with observations one faces 
a difficult problem — the overwhelming majority of these 
observations do not properly isolate the torus IR emission. 
Starburst emission is increasingly recognized as an important 
component of the IR flux me asured in many, perhaps most, 
AGN (e.g jNetzer et all2007l) . In addition to this well known 
contamination, even IR from the immediate vicinity of the 
AGN may not always originate exclusively from the torus, 
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furt her complica t ing m odeling efforts. A case in point is 
the iMason et all rf2006h modeling of NGC1068. All flux 
measurements with apertures < 0.5" are in good agreement 
with the model results, but the flux collected with larger aper- 
tures greatly exceeds the model predictions at wavelengths 
longer than ~4fim. This discrepancy can be attributed to IR 
emission from nearby dust outside the torus. Mason et al 
show that the torus contributes less than 30% of the 10 /im 
flux collected with apertures > 1" and that the bulk of the 
large-aperture flux comes at these wavelengths from dust in 
the ionization cones; while less bright than the torus dust, 
it occ upies a much larger volume (see also iPoncelet et al.l 
120071) . On the other hand, the torus dominates the emission 
at short wavelengths; at 2 fj,m, more than 80% of the flux 
measured with apertures > 1" comes from the torus eve n 
though its image size is less than 0.04" (Weigelt et al. 2004). 

These difficulties highlight a problem that afflicts all IR 
studies of AGN. The torus emission can be expected to dom- 
inate the AGN observed flux at near IR because such emis- 
sion requires hot dust that exists only close to the center. But 
longer wavelengths originate from cooler dust, and the torus 
contribution can be overwhelmed by the surrounding regions. 
Unfortunately, there are not too many sources like NGC1068. 
No other AGN has been observed as extensively and almost 
no other observations have the angular resolution necessary 
to identify the torus component, making it impossible to de- 
termine in any given source which are the wavelengths dom- 
inated by torus emission. There are no easy solutions to this 
problem. One possible workaround is to forgo fitting of the 
spectral energy distribution (SED) in individual sources and 
examine instead the observations of many sources to iden- 
tify characteristics that can be attributed to the torus signa- 
ture. One example for the re moval of the starburst compo- 
nent is the lNetzer et al.1 (12007) composite SED analysis of the 
Spitzer observations of PG quasars. Netzer et al identify two 
sub-groups of "weak FIR" and "strong FIR" QSOs and a third 
group of far-IR (FIR) non-detections. Assuming a starburst 
origin for the far-IR, they subtract a starburst template from 
the mean SED of each group. The residual SEDs are remark- 
ably similar for all three groups, and thus can be reasonably 
attributed to the intrinsic AGN contribution, in spite of the 
many uncertainties. However, while presumably intrinsic to 
the AGN, it is not clear what fraction of this emission origi- 
nates from the torus as opposed to the ionization cones. An 
example of samp le analysis t hat m ay have identified the torus 
component is the lHao et al.l d2007l) compilation of Spitzer IR 
observations. In spite of the large aperture of these measure- 
ments, Seyfert 1 and 2 galaxies show a markedly different be- 
havior for the 10/im feature, both in their mean IR SEDs and 
in their distributions of feature strength. Furthermore, Ultra- 
luminous IR Galaxies (ULIRG) that are not associated with 
AGN show yet another, entirely different behavior, indicating 
that the observed mean behavior of Seyfert galaxies is intrin- 
sic to the AGN. Accepting the framework of the unification 
scheme, the differences Hao et al find between the appear- 
ances of Seyfert 1 and 2 can be reasonably attributed to the 
torus contribution; the ionization cones dust is optically thin, 
therefore its IR emission is isotropic and cannot generate the 
observed differences between types 1 and 2. 

Here we invoke both approaches in comparing our model 
predictions with observations. We start by assembling dusty 
clouds into complete models of the torus, as described in $2] 
Our model predictions for torus emission and the implications 
to IR observations are presented in §§|3rll] while in $6]we dis- 




FlG. 1. — Model geometry: Dusty clouds, each with an optical depth tv 
at visual, occupy a toroidal volume from inner radius Kj, determined by dust 
sublimation (eq.fTJ, to outer radius R„ = YR L \. The radial distribution is a 
power law r~ q , the total number of clouds along a radial equatorial ray is 
Afo- Various angular distributions, characterized by a width parameter cr, 
were considered. The angular distribution has a sharp-edge on the left and a 
smooth boundary (e.g., a Gaussian) on the right. 

cuss aspects of dumpiness that are unrelated to the IR emis- 
sion, such as the torus mass, unification statistics, etc. In $7] 
we conclude with a summary and discussion. 

2. MODEL OF A CLUMPY TORUS 

Consider an AGN with bolometric luminosity L surrounded 
by a toroidal distribution of dusty clouds (Fig. [T]). The 
"naked" AGN flux at distance D is Fagn = L/AttD 2 at any 
direction, but because of absorption and re-emission by the 
torus clouds the actual flux distribution is anisotropic, with 
the level of anisotropy strongly dependent on wavelength. 
The grain mix has standard interstellar properties (see paper 
I, §3.1.1 for details), the optical depth of each cloud is Ty at 
visual. 

2.1. Dust Sublimation 

The distribution inner radius Ra is set by dust sublimation 
at temperature T su b. From §3. 1 .2 in part I, 

^^0.4 —^ -A — pc. (1) 

\10 43 ergs V V Tsub / 

iBarvainisI (119871) derived an almost identical relation for R4. 
His eq. 5 has the same normalization and only a slight dif- 
ference in the power of T su \, (2.8 instead of 2.6); this dif- 
ference reflects the more detailed radiative transfer calcula- 
tions we perform. Here the distance is determined from 
the temperature on the illuminated face of an optically thick 
cloud of composite dust representing the grain mixture. The 
sharp boundary we employ is an approximation. In reality, 
the transition between the dusty and dust-free environments 
is gradual because individual components of the mix subli- 
mate at slightly differen t radii, with the largest g rains surviv- 
ing closest to the AGN (ISchartmann et al.ll2005l). From near- 
IR reverberation me asurements, iMinezaki et alj d2004l) and 
Sugan uma et alj d2006l) find that the inner radius of the dusty 
region is indeed proportional to L 1 / 2 , but the time lags they 
report are ~2-3 times shorter than predicted by eq. [T] While 
this equation gives the smallest radius at which the dust ab- 
sorption coefficient reflects the full grain mixture, the largest 
grains survive to closer radii, where they are presumably de- 
tected by the reverberation measurements. 

2.2. The Cloud Distribution 

The torus extends radially out to R Q = YR^, with Y a free 
parameter. The total number of clouds, on average, along any 
radial equatorial ray is specified by the parameter No. We 
studied various forms for the variation of A/t(/3), the total 
number of clouds along rays at angle from the equator. Fig- 
ureQ]shows on the left a sharp-edge uniform distribution with 
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Viewing Angle, i (deg) 

FIG. 2. — Behavior of the probability for photon escape along a path con- 
taining Af clouds (eq. 4 part I). Upper panel: Wavelength variation of P e sc for 
the indicated Af when the single cloud optical depth is ry at visual. Lower 
panel: The probability for an AGN photon to escape through the torus in 
direction i from the pole when each cloud is optically thick; this is also the 
probability for unobscured view of the AGN at viewing angle i. The to- 
tal number of clouds varies according to A/t(/3) = A/o exp(— 1 /a 1 ), where 
j3= bTT— lis angle from the equatorial plane, with a = 45° and A/o as marked. 

A/t(/3) = A/o within the angular width \(3\ <a. In a Gaussian 
distribution, A/ T (/3) = Af exp(-f3 2 /a 2 ). 

The emission from the clumpy torus is found by integra- 
tion along paths through the cloud distribution (eq. 5, part I). 
Some of the computation technicalities are described in the 
Appendix. The calculation requires the single cloud source 
function S Ci \, derived in part I, and the number of clouds per 
unit length, Nc(r,(3), as a function of (3 and radial distance 
r. For this distribution we assume a separable function with 
power law radial behavior r~ q so that 



N c (r,(3) = C — — x — 
fid V r 



(2) 



where C = (J; dy/y^T 1 is a dimensionless constant (for 
a given Y and a), ensuring the normalization Nr(ff) = 
jN c (r,f3)dr. 

The observed torus radiation is affected not only by the 
emission from individual clouds but also by the probability 
that emitted photons escape through the rest of the path. The 
escape probability, P e sc, is given in eq. 4 part I. For an overall 
number of clouds M along a path, P esc ~ exp(-A/V\) at wave- 
lengths in which the optical depth of a single cloud obeys 
r\ < 1 , and P esc ~ exp(-A0 when t\ > 1 . Many of the de- 



tailed results presented below can be readily understood from 
the dependence of P esc on wavelength and on torus viewing 
angle, shown in figure [2] 

2.3. Scaling 

Be cause of general scalin g properties of radiatively heated 
dust dlvezic & Elitzur|[l997l) . the only effect of the overall lu- 
minosity is in setting up the bolometric flux F agn and the dust 
sublimation radius /?d (eq. [TJi. For a given torus model, the 
distributions of dust temperature and of brightness are unique 
functions of the scaled radial distance r/R&: two sources with 
the same cloud properties but different luminosities will have 
the same distributions in terms of r/R^, only the more lumi- 
nous one will have its brightness spread over a larger area 
because of its larger R& (this point is explained further in the 
Appendix). Denoting the torus flux by F\, the spectral shape 
F\/Fagn is independent of L. The dependence of the torus 
SED on the spectral shape of the AGN input radiation is lim- 
ited to scattering wavelengths, disappearing altogether at A > 
2-3 /mi. There is a similarly weak dependence on r su b. The 
output spectrum depends primarily on -ry and the cloud distri- 
bution. Although the luminosity does not affect the radiative 
transfer, it is entirely possible for torus properties to be corre- 
lated with L for some other reasons (e.g. a, as in the receding 
torus model). 

2.4. The AGN Contribution 

In most figures we show only the contribution of the torus 
emission. However, since the medium is clumpy, there is al- 
ways a finite probability for an unobscured view of the AGN, 
irrespective of the viewing angle. Because of the probabilistic 
nature of the problem it is only possible to display the emerg- 
ing spectral shape with or without the AGN contribution and 
the probability for each case (see bottom panel of figure[2]i. 

3. MODEL SPECTRA 

We proceed now with the model results. In all calculations 
the AGN input radiation follows the "standard" spectrum de- 
scribed in §3.1.1, part I. 

3.1. Geometrical Shape 

Figure[3]shows model results for sharp-edged and Gaussian 
angular distributions. The sharp-edge geometry produces a 
bimodal distribution of spectral shapes, with little dependence 
on viewing angle other than the abrupt change that occurs be- 
tween the torus opening and the obscured region. In contrast, 
the Gaussian distribution produces a larger variety in model 
spectral shapes, with a smooth, continuous dependence on i. 
We investigated a larger family of angular distributions of the 
form A/t(/3) = A/oexp(— \(3/a\ m ), with m a free parameter. In 
this family, m = 2 is the Gaussian and as m increases, the 
transition region around (3 = a becomes steeper. Generally, 
"softer" distributions with m < 10 show behavior similar to 
the Gaussian while those with larger m produce results simi- 
lar to the sharp-edge geometry. 

The SED dichotomy produced by sharp boundar ies con- 
flicts with observations. Alon so-Herrero et"afl (120031) studied 
the 0.4-16 /im nuclear emission from a complete sample of 58 
Seyfert galaxies, selected fro m the CfA sample. In a compar- 
ison with theoretical models, lAlonso-Herrero et alJ point out 
that a common prediction of all smooth-density models is a di- 
chotomy of SED between type 1 and 2, similar to the one dis- 
played in the upper panel of fig. [3] and that such a dichotomy 
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FIG. 3. — Model spectra for a torus of clouds, each with optical depth 
Ty = 60. Radial distribution with q = 1 out to Y = 30, with JVo = 5 clouds 
along radial equatorial rays (see eq. ff). The angular distribution is sharp- 
edged in the top panel, Gaussian in the bottom one (cf. fig. [T}; both have a 
width parameter cr = 45°. Different curves show viewing angles that vary 
in 10° steps from pole-on (7 = 0°) to edge-on (i = 90°). Fluxes scaled with 
Fagn=L/4itD 2 . 

is not observed in their sample; the dichotomy is present even 
in model geometries with soft edges because the exp(-r) at- 
tenuation factor varies rapidly, resulting in a sharp transition 
around r ~ 1 between dusty and dust-free viewing. As is 
evident from the lower panel of fig. [3] this SED dichotomy 
problem is solved by soft-edge clumpy tori. Therefore, in the 
following we consider only Gaussian angular distributions. 

3.2. Observations and Model Parameters 

As discussed in the Introduction, torus IR observations are 
hampered by uncertainties that are partially alleviated by con- 
sidering composite spectta. Figure |4] shows compilations of 
type 1 and type 2 data and some r epresentative model s, up- 
dating a similar figure presented in Nen kova et al.l d2002l) . The 
type 1 data a dditionally includ e th e recent Spitz e r com posite 
spectta from lHao et ail (120071) and lNetzer etail (120071) . The 
close agreement between these two SEDs in their common 
specttal region, A = 5-3 8 /im, indicates that they may have 
captured the torus emission in outline, if not in details. The 
upturn around 60/im in the Netzer et al spectrum likely re- 
flects the transition to starburst dominance. To ensure the 
smallest possible apertures in type 2 sources, the data for in- 
dividual objects are mostly limited to ground based and HST 
observations. The data in both panels of this figure display 
the general characteristics that have to be reproduced by the 
same models in pole-on and edge-on viewing. The updated 
models plotted with the data differ from the original ones in 
iNenkova et al.l (T2002) in three significant ways: (1) the optical 
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FIG. 4. — Observations of type 1 and type 2 sources compare d with clumpy 
torus model spectra. Th e type 1 composi te da ta are from I Sanders et all 
fl989l) . Elvis et al.1 (19941) . IHao et alj f200l a nd INetzer et alj 120071) . The 
type 2 data are from the following sources: a) IMason et alj 120061); b) var- 
ious observations wi t h aper ture < 0.5" lis ted in Mason et al. 12006); c) 
lAlonso-Herrero et al.1 <2003l) ; d) Prieto et afj 120041) . In the model calcula- 
tions, plotted with broken lines, each cloud has optical depth Ty = 30. Other 
parameters are cr = 30°, q — 0-3, as marked, Y = 30 and A/o = 5. The angular 
distribution in this and all subsequent figures is Gaussian. The models in the 
upper panel are for pole-on viewing ((' = 0° ), in the bottom panel for edge-on 
(i = 90°). 

properties of the silicate component o f the dust are taken from 
the tabulatio n for "cool" sil i cates in lOssenkopf et alj (119921) 
instead of the lDraine & Led (|T984) dust; (2) the clouds angu- 
lar distribution is Gaussian rather than sharp-edged; and (3) 
the torus radial thickness Y is 30 instead of 100. As is evident 
from the figure, the model spectra are generally in reasonable 
agreement with the data. 

We produced a large number of models for various param- 
eter sets 4 , and we now present model results and discuss their 
observational implications. The models are characterized by 
free parameters that describe individual clouds (tv), control 
the total number of clouds (A/o) and specify the geomettical 
properties of their angular and radial spatial distributions (cr, 
q and Y). Note that, except for Ao, smooth-density models 
also require all of these parameters to describe the dust dis- 
tribution. In the following, the parameters are varied one at 
a time, and from comparison with observations we attempt to 
identify the likely range of each of them. The effect of the ra- 
dial thickness parameter Y is described separately in the next 
section (®, devoted to a discussion of the torus size. 

3.3. Single Cloud Optical Depth 

4 Tabulations of all the models discussed here as well as many additional 
cases are available at |http : / /www .pa ■ uky . edu/ clumpy /| 
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FIG. 5. — Dependence of the torus SED on the single cloud optical depth 
tv . Other parameters are a — 45° , A/o = 5 and Y = 30. Radial power law with 
q — 1 in the left panels, q = 2 in the right ones. Pole-on viewing in the top 
panels, edge-on at the bottom. 

Figure [5] shows the effect of varying the optical depth of in- 
dividual clouds from 10 to 200. The SED hardly varies when 
Ty increases beyond ~ 100, reflecting the similar behavior for 
emission from a single cloud (see Fig. 11, Part I). The fig- 
ure shows the torus emission for both pole-on and edge-on 
viewing. Smooth-densit y models (e.g.. iPier & Kjoliklll992l : 
iGranato & Danesel [T994) consistently produce the 10 /im sili- 
cate feature in emission and absorption, respectively, for polar 
and equatorial viewing. As the figure shows, in a clumpy dis- 
tribution the feature displays a more complex pattern, unlike 
anything produced in smooth-density models. At i = 0°, the 
feature appears in emission as long as Ty < 20. When the 
optical depth increases further, the feature disappears and the 
SED is essentially featureless across the 10/im region. How- 
ever, the feature reappears in weak emission when Ty > 100. 
At i = 90°, a weak, broad emission feature is evident when 
Ty = 10. When Ty > 20, the spectra display a clear absorp- 
tion feature; although similar to that of smooth-density mod- 
els, the feature is never deep, reflecting the shallow absorption 
displayed by a single cloud (see §4.5 part I). A most peculiar 
result is the reversal from absorption to an emission feature, 
which emerges when Ty increases beyond ~ 100. 

The complex behavior of the 10 /im feature arises from a 
rather intricate interplay between the emission spectrum of a 
single cloud and the collective effect of the entire cloud en- 
semble. The different patterns can be understood in terms 
of the competition between emission and absorption along a 
given path, taking account of the flattening of the escape fac- 
tor P esc across the 10/im feature when Ty is increasing (fig. 
13. The behavior of th e 10 /im feature is studied separately at 
greater depth in {35.1l below. 

3.4. Number of Clouds 

Figure [6] shows model spectra of torus emission when A/o, 
the average of the total number of clouds along radial equato- 
rial rays, varies from 2 to 15. The models produce broad IR 
emission in the A ~ 1-100 /im range. Values of A/"o larger than 
15 produce a very narrow IR bump peaking beyond 60 /im. 
Such SED's have not been observed thus far, therefore A/o is 
likely no larger than ^10-15 at most. 

As is evident from figure [6] when A/o increases, the emis- 
sion in the near- to mid-IR region steepens considerably for 
viewing close to equatori al. Composite IR SED of Seyfert 
2 galaxies constructed by ISilva et al.l ([2004) show only mild 
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FIG. 6. — Dependence of the torus SED on the number of clouds along a 
radial equatorial ray. Each cloud has Ty = 60. Angular width a = 45°, power 
law radial distribution with q = 1 (left panels) and 2 (right), extending to Y 
= 30. Different curves in each panel show viewing angles that vary from 0° 
(top curve) to 90° (bottom) in 10° steps. 
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FIG. 7. — Model spectra as in figure |6]for q = 2, only with the AGN 
contribution added. For each set of parameters, the probability that the AGN 
emission will actually be observed can be read from the plots of P e sc in the 
lower panel of figure [2] The break in the SED at A = 1 is an artefact of 
our parametrization of the input spectrum (see §3.1.1, part I). 



dependence on X-ray absorbing column density as long as 
Nh < 10 24 cm -2 , with considerable steepening when Nn is 
in the range 10 24 -10 25 cm" 2 . This behavior is similar to the 
A/()-dependence displayed in fig. [6] thus Compton thick X-ray 
absorbing columns might be correlated with a larger No. For 
pole-on viewing, the 10 /im feature appears in weak emission 
when No < 5. As No increases, the emission switches to ab- 
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FIG. 8. — Dependence of the torus model spectra on the width parameter 
a of the Gaussian angular distribution. Each cloud has -ry = 60. The cloud 
radial distribution is a power law with A/"o = 5 and q = 2 extending to Y = 30. 
Viewing angles vary from 0° to 90° in 10° steps. 

sorption that deepens with JVq. Moving away from the axis, 
the feature displays weak emission when A/o = 2 but appears 
in absorption in all other cases. 

In contrast with the smooth-density case, clumpy models 
always display some emission at A < 1 /im that arises from 
scattering of the AGN radiation toward the observer by clouds 
on the torus far side. Some fraction of this radiation will al- 
ways get through the torus near side. The probability for that 
is controlled purely by the number of clouds since individ- 
ual clouds are always optically thick at UV and optical wave- 
lengths. Varying the number of clouds produces two compet- 
ing effects, most clearly visible in figure|6]from the behavior 
of the q = 1 models at i = 0°. Increasing the number of clouds 
from A/o = 2 to 5 raises the level of the radiation that gets 
through because there are more scattering clouds. With fur- 
ther increase in Ao, obscuration by intervening clouds takes 
over and the emerging intensity decreases. It is hard to as- 
sess the observational significance of this aspect of the results. 
Our models include a single type of clouds and no intercloud 
medium. Such a medium with an optical depth -ry of only a 
few would attenuate all wavelengths shorter than ~1 /im in 
the model spectra without significantly affecting the infrared. 
We plan a detailed study of these effects in future work. 

Like most models presented here, fig. [6] shows only the 
torus emission, corresponding to type 2 SED. Our model pre- 
dictions for type 1 SED can always be obtained by adding the 
AGN direct radiation. However, unlike the smooth-density 
case, the probability for a clear view of the AGN depends not 
only on the viewing angle but also on the number of clouds 
(see jj2.4| i. Figure Q displays again the q = 2 models shown 
in fig. [6] but this time the AGN contribution is added in. The 
probability that this would be the SED actually detected in a 
given source is given by the corresponding P esc , shown in the 
lower panel of figure [2] When visible, the AGN dominates 
the emission at A < 3/im. The transition from AGN to torus 
domination of the SED is an important issue that requires de- 
tailed observations of type 1 sources in the near- and mid-IR 
regions. 

3.5. The Torus Angular Width 

The effect of the angular distribution width is shown in fig- 
ure [8] which displays results for a few representative a. The 
spectral shapes of models with a = 15° are in general agree- 
ment with observed SED but the dependence on viewing an- 
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FIG. 9. — Dependence of the torus model spectra on the power q of the 
radial density distribution, which extends to Y = 10 in the top panels and Y = 
30 in the bottom ones. A/"o = 5 clouds with Ty = 60 each. Angular width cr = 
45°. Viewing angles from 0° to 90° in 10° steps. The emission anisotropy 
decreases when q increases. 

gle displays a bi-modal distribution that conflicts with obser- 
vations of Seyfert galaxies Q3.lt . Values in the range 30°- 
50° produce similar spectral shapes, all in general agreement 
with observations. The a = 30° models provide the best match 
to the behavior of the lOfim fea ture in the average spectra of 
Seyfert 1 and 2 galaxies (see §15. 11 1. Estimates of the torus 
angular width based on statistics of Seyfert galaxies that take 
proper account of dumpiness give a ~ 30° (see ^6.3t . At a 
= 60°, the lOfim feature appears in pronounced absorption at 
all viewing angles. Increasing the width parameter further all 
the way to a = 85° has little effect on the SED, except that the 
dependence on viewing angle decreases, as is expected from 
the approach to spherical symmetry. 

3.6. Radial Profile and IR Emission Anisotropy 

Figure [9] shows the SED when q, the index of the power- 
law radial distribution, varies in the range 0-3 for two values 
of the radial thickness Y. Since A/o is kept fixed, varying q 
changes only the placement of clouds between the inner, hot- 
ter parts and the outer, cooler regions, shifting the emission 
between near- and far-IR. Steep radial distributions (q = 2, 3) 
produce nearly identical spectral shapes for Y = 10 and 30 be- 
cause the clouds are concentrated near the inner boundary in 
these cases. 

The variation of SED with viewing angle displayed 
in fig. |9] is much smaller t han in smooth-density models 
jPier &^rolikl 1 19921 119931; lEfstathiou & Rowan-Robinsonl 
1 1 9951: iGranato & DaneseT Il994t iGranato et alj Il997t 
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FIG. 10. — Variation of the torus flux with viewing angle at different wave- 
lengths, as marked. JVo = 5 clouds with Ty = 60 each in a radial distribution 
with q = 1 (top panel) and 2 (bottom) extending to Y = 30. Angular width a = 
45°. Since the emission is normalized to the AGN flux, the plotted quantity 
provides the bolometric correction for each displayed wavelength. 



iDullemond & van Bemmel 120051: ISchartmann et all 120051) . 
For example, our models produce at 1 /mi an edge-to-pole 
flux ratio of abou t 5 or less for cloud optical depth of Ty = 60 
. In Fig. 1 from iGranato et all the corresponding flux ratio 
is several hundred for Ay = 30, and off-scale for Ay = 100. 
The degree of isotropy is especially high when the torus is 
small, but even in the case of Y = 30 the emission becomes 
nearly isotropic as q increases (the radial distribution gets 
steeper). This point is further illustrated in Figure [TO) which 
shows the viewing-angle variation of the observed flux at 
different wavelengths. Since the torus flux is normalized 
to the AGN overall flux, the quantity plotted in this figure 
provides the bolometric correction for each of the displayed 
wavelengths. The anisotropy decreases with wavelength, 
practically disappearing beyond the mid-IR — at A = 12 
/im the variation with viewing-angle is within a factor of 
^2 for both q = 1 and 2. Another indicator of the emission 
anisotropy is the variation of the torus bolometric flux F tor = 
J F\d\ with viewing angle, shown in figure[TT]for q = 2. The 
variation increases with the number of clouds but remains 
less than a factor of 3 even when Afo = 15, the likely upper 
limit. 

It is important to note that at every viewing angle, the AGN 
obscuration is identical in all the models displayed in fig. |9l 
obscuration depends only on the total number of clouds along 
radial rays, which is the same in all cases. Indeed, for each of 
these models the probability for direct view of the AGN as a 
function of ; is shown by the Afo = 5 curve in the bottom panel 
of figure |2] Along this curve, P esc varies by two orders of 
magnitudes between polar and equatorial viewing. That is, a 
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FIG. 1 1. — Anisotropy indicators for the torus bolometric flux F tm . All 
models have q = 2 and Y = 30. Top: Variation of F t0I with viewing angle 
when Ty = 60 for various values of JVo (cloud number) and <r (torus angular 
width), as indicated; note the changing vertical scale. Bottom: Ratio of the 
torus bolometric fluxes along the axis and equator as a function of single 
cloud optical depth. The torus angular width is cr = 45° . 

clumpy torus can produce extremely anisotropic obscuration 
of the AGN together with nearly isotropic mid-IR emission. 

Recent observations seem to indicate that this is indeed 
the required behavior. Ground based observations of AGN 
nuclear emission at 10/i,m show it to be well correlated 
with the hard X-ray luminosity, and both type 1 and type 
2 sources follow the sam e correlation dAlonso-Herrero et alj 
l200UlKrabbe et alj|2001l) . IWhvsong & Antonuccil d2004 re- 
port that the mean values for the 12 /im/60 /im flux ratios of 
Seyfert 1 and 2 galaxies differ by only ~30%. Since the 60 
im\ emission is optically thin and thus essentially isotropic 
(cf fig. [T0l). t his result indicate s that the variation of the 12 /im 
flux is small. iLutz et all (H004) compared the 6 /im ISO fluxes 
of Seyfert 1 and Seyfert 2 galaxies normalized to their intrin- 
sic hard X-ray fluxes. They conclude that the distributions 
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FIG. 12. — Dependence of the SED of a clumpy torus on the radial thickness 
Y = R /Rd, as marked. Radial distribution with q = 1 (left panels) and 2 
(right). All models have JVo = 5 clouds with ry = 60 each, and cr = 45°. Pole- 
on viewing in the top panels, edge-on at the bottom. Note that the curves in 
the left-bottom panel have a similar shape at A < 15 /im and would nearly 
overlap if normalized to a common wavelength in that range instead of Fagn- 

of the two populations are essentially identical within the ob- 
servational errors, and note the conflict with the anisotropy 
predicted by smooth-dens ity torus models. The iLutz et all 
finding was confirmed by iHorst et all ([2006), who used the 
same approach with grou nd based, and thu s better resolution, 
observations at 12 /im. [Buchanan et al.l (120061) conducted 
Spitzer observations of 87 Seyfert galaxies in the A = 5-35 
/im range and normalized the IR fluxes with the optically thin 
radi o emission. Although at 6 /im they find a larger variation 
than lLutz et al.L they also find that the emission from Seyfert 
1 and 2 galaxies are within factor 2 of each other for all A > 
15 /im, and note the discrepancy with smooth-density mod- 
els. Finall y, the aver a g e spectra of Seyfert 1 and 2 galaxies 
derived by [Hao et all d2007l) from Spitzer observations have 
nearly identical shapes, except for the 10/im silicate region. 

The moderate level of anisotropy found in the observations 
suggests that, if the torus radial thickness is > 20 then the 
steeper radial profile q = 2 might be more appropriate than 
<7=1. It may be noted that the clumpy torus models in 
iMason et all (120061) . which utilized Y = 100, yielded the best 
fits to the observations of NGC 1068 with q = 2. 

4. TORUS SIZE 

The fraction of the sky obscured by the torus determines the 
relative numbers of type 1 and 2 sources, and the statistics of 
Seyfert galaxies show that the h eight and radius of obscuring 
dusty torus obey H/R ~ 1 (see 96.31 . Since obscuration does 
not depend separately on either H or R, only on their ratio, 
neither quantity is determined individually. An actual size 
can only be determined from the torus emission. 

4.1. SED Analysis 

In the absence of high-resolution IR observations, early es- 
timates of the torus size came from theoretical analysis of the 
SED. For a given dust sublimation temperature, the torus in- 
ner radius R^ is determined from the AGN luminosity (eq. 
[T|). The dust temperature distribution, and with it all model 
results, depends only on r/R^. Therefore, the only size pa- 
rameter that can be determ ined from SEP mod eling is the 
radial thickness Y = R /Rd- iPier & Krolikl (|T992) performed 
the first detailed calculations with a uniform de nsity torus and 
found Y ~ 5-10. However, in subsequent work lPier & Krolikl 



( 1 1993b speculated that this compact structure might be em- 
bedded in a much larger, and more diffuse t orus, extending 
typically to - 30-100 pc. lGranato & Danesdd 19941) extended 
the smooth-density calculations to more elaborate toroidal ge- 
ometries. They conclude that "The broadness of the IR con- 
tinuum of Seyfert 1 nuclei requires an almost homogeneous 
dust distribution extending at least to a few hundred pc (R /Rd 

> 300 or R Q > 300L^ 2 pc)" and that "broad (R Q ~ lOOOL^ 2 
pc) tori" would be "fully consistent with available broad-band 
data and high-resolution IR spectra of Seyfert 1 and 2 nu- 
clei". Althoug h subsequent modeling produced somewhat 
smaller sizes dGranato et al.lll997t Fritz et al.ll2006l) . the orig- 
inal requirements of uniform density and large dimensions 
directly reflect the large amounts of cool dust necessary for 
producing the torus IR emission. This requirement arises be- 
cause in smooth density distributions, the dust temperature 
is uniquely related to distance from the AGN. While this 
statement is strictly correct only for single-size dust grains, 
even when dust size distribution in invoked in smooth dust 
models, the observation s still favor clumpy dust distribution 
(ISchartmann et alJl2005h . 

The one-to-one correspondence between distance and tem- 
perature does not hold in clumpy media, where different dust 
temperatures coexist at the same distance and where the same 
temperature can be found at different distances (see part I, 
§3.1.2 and §4.2). For example, in the model discussed in 
Fig. 7 from part I, the dus t temperature at Y =10 ranges from 
150 K to 600 K, while ISchartmann et alJ ((2005) find using 
smooth models with realistic dust size distribution that the 
temperature range at Y = 10 is 250-300 K (i.e. a ratio of 1.2 
vs. 4 for clumpy dust). In contrast with smooth density dis- 
tributions, a clumpy torus contains cool dust on the dark sides 
of clouds much closer to the heating source and thus can emit 
IR efficiently from its inner regions. 

Figure [T2] shows our model results for the SED of clumpy 
tori with various radial thicknesses for q = 1 and 2. In spite of 
the factor 40 variation in torus thickness, the SED is quite 
similar for all the q = 2 models. The reason is simple — 
irrespective of the torus size, at least 80% of all clouds are 
located within r < 5R& in this case. The models with q = 1 
display discernible variations, but these variations are mostly 
confined to A > 15/im. The large differences apparent at 
shorter wavelengths in the edge-on viewing do not reflect in- 
trinsic variation of the SED, only the scaling with the bolo- 
metric flux. If these curves were scaled instead to the same 
value at, say, 2/im, they would all overlap up to A ~ 15/im, 
similar to the pole-on viewing. These differences can be fur- 
ther understood with the aid of figure [13] which shows the 
fraction of the overall flux contained within circular apertures 
of increasing size. This fraction, as well as the brightness 
distribution, is a function of 6 /6a (= r/R^), where 8 is angu- 
lar displacement from the center and f?a = Ra/D. The figure 
shows that at wavelengths shorter than 5 /im, almost all the 
flux is originating from inside 36 a irrespective of the value of 
Y. Therefore, such wavelengths cannot determine the torus 
size. In the q = 2 case, even longer wavelengths cannot dis- 
tinguish between the different sizes because 80% of the flux 
always originates from the inner 10 6*d- In the q = 1 case the 
portions beyond 106 a contribute significantly to the flux of a 
larger torus, but only at wavelengths longer than ^12-15/im. 

These results show that the model SED do display apprecia- 
ble differences among tori of different sizes when q = 1 and 
that determining the torus size in this case would require mea- 
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FIG. 13. — Fraction of the torus flux enclosed within a circle with angular radius 8 centered on the AGN. Wavelengths, in (im, as labeled (some labels omitted 
for clarity). 9$ is the angle equivalent of the torus inner radius R$ (eq.[T) at the observer's location; for reference, S ( j= 0"02 for i?j = 1 pc at 10 Mpc. All models 
have A/o = 5, 7y = 60, a — 45° and ^ = 1 or q = 2, as marked. In each case, pole-on viewing in the left panels, edge-on at right. Torus sizes, from top to bottom, 
are Y = 10, 30 and 100. 



surements of its flux at A > 15 /im. Because of the large beam 
sizes at such long wavelengths, current observations generally 
cannot distinguish between the contributions of the torus and 
its surroundings to the overall flux. For example, as noted 
in the Introduction, in NGC 1068 the torus contributes less 
than 30% to the 1 /im flux measured with apertures > l" 
( Maso n et alj2006l) . The q = 2 density profile, which might be 
the more common radial distribution Q3.6b . does not gener- 
ate any discernible distinctions among spectral shapes. There- 
fore, the radial size of a clumpy torus cannot be constrained 
by SED measurements. Only high-resolution observations can 
determine the torus size — the SED does not have the neces- 
sary discriminative power. 

4.2. Brightness Profiles 

Because of the discrete nature of a clumpy medium, differ- 
ent lines of sight will generally produce a different brightness 
even when crossing similar regions. Our formalism provides 
the statistical average along such rays; fluctuations around this 
average can be large (see §2.1, part I). 

Figure [14] shows model intensity profiles for pole-on view- 
ing of a torus with Y = 30. The curves show only the torus 
emission, starting at its inner edge 9 = 9&. The AGN emission, 
which is not shown, would produce a narrow spike between 
< 9 < 9 AGN , where 9 AGN /9 d ~ (7 Wr b , A GN) and T t ,.agn i s 
the AGN brightness temperature (e.g.,Hyezic & El itzurl 1997b . 
Therefore, 9 A gn / 9 d <C 1 under all circumstances. The set 
of displayed wavelengths extends from the /^-band, where 
most of the current imaging observations are performed, to 
588 /im, one of the wavelengths that will become available 
for high-resolution imaging when ALMA is fully operational. 
The torus intensity is highest on or close to its inner edge. 
The brightness is highest around 12 /im, as is evident from 



the figure bottom panel. For both radial density profiles used 
in this figure, the brightness declines to half its peak value 
within 9 < 59 d at all displayed wavelengths. At 12/im and 
shorter wavelengths, the brightness declines to 1% of peak 
value within 9 < 10 6*d- Evidently, observations attempting to 
probe the torus structure must combine high resolution with a 
large dynamic range. 

Near-IR wavelengths provide little information about the 
torus structure and size. As is evident from the figure, at 2 /im 
it would be difficult to distinguish between q = 1 and 2 radial 
density distributions even with high-resolution observations. 
The steep brightness decline at these wavelengths also makes 
it practically impossible to determine the torus full size. Since 
the brightness falls under 1 % of its peak value at 9 = 4 9 d for 
either density profile, determining whether the torus ends at 
that point or continues to larger radii would be a difficult task. 
A Y = 10 torus is indistinguishable from the inner 106>d of a 
torus as large as Y = 100, as is evident also from figure [13] 
As the wavelength increases, the brightness fall-off becomes 
less steep. VLTI interferometry has angular resolution of or- 
der Q'.'Ql at 12 /im, but it would still be difficult to distinguish 
between the two displayed radial density profiles even in sys- 
tems where #d is of a similar order of magnitude. The two 
density profiles produce distinctly different brightness profiles 
at 70/im, but there are no instruments with the required angu- 
lar resolution at that wavelength. In the foreseeable future, 
ALMA seems to be the only facility with a realistic chance 
to determine through 588 /im observations the radial cloud 
distribution and whether a torus does extend beyond Y = 30. 

4.3. Observations 

With the advent of high-resolution IR observations, di- 
rect imaging is now available for some AGN tori, and up- 
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FIG. 14. — Surface brightness for torus models with -ry = 60, A<o = 5, a 
= 45°, Y = 30 and q = 1 and 2, as indicated. The top panel shows the radial 
variation of intensity with angular displacement from the center for pole-on 
viewing and a set of wavelengths as marked. The AGN emission, which is not 
shown, corresponds to a narrow spike at 9/9^ <C 1 (see text). Each intensity 
profile is normalized to its brightness level at 9 = 8i, shown in the bottom 
panel together with the corresponding brightness temperature. 

per limits have been set on the dimensions of the nuclear IR 
source in others. Interferometric observations at 8-13 /Ltm 
with the VLTI have resolved by now the nuclear region in 
three AGN: NGC 1068, Circinus and Cen A. The thermal 
emission in all t hree cases is rather compact. In NGC 1068 
Jaffe et al. (2004) find that the emission extends to R = 1.7 pc. 
Poncelet et al.1 (2006) reanalyzed the same data with slightly 
different assumptions and obtained a similar result, R = 2.7 
pc. The AGN bolometric lu minosity is ~ 2xl0 45 ergs" 1 in 
this case dMason et al.ll2006T) . so that /?d is —0.6 pc and the 
toru s mid-IR emission is confined within ~3— 5R&. In Circi- 
nus, |TOsttametal](|2007) find that the torus emission extends 
to R = 1 pc. The AGN bolometric luminosity is ~ 8xl0 43 
erg s -1 (lOliva et al.l fT999). so R^ — OA pc and the outer radius 
of the mid-IR emission is — 10/?d- The nature of the mid- 
IR emission from the Cen A nucleus is somewhat involved — 



iMeisenheimer et ail d2007l) conclude that it contains an un- 
resolved synchrotron core and thermal emission within a ra- 
dius of —0.3 pc . Since the AGN bolometric luminosity is ~ 
1 x 10 43 erg s" 1 dWhvsong & Antonu cci 2004), R d is -0.04 pc 
and the torus emission does not exceed — 8/?d in this source. 
One other c ase of resolved mid -IR emission involves NGC 
7469, where ISoifer etail d200l find a 12.5 fim compact nu- 
clear structure contained within R < 13 pc. Unfortunately, 
NGC 7469 is a clear case where the IR signature is dominated 
by the starburst compone nt even though the AGN dominates 
the optical classification ( Weedman et al. 2005), therefore the 
resolved compact structure c annot be identified with the torus 
(see also lDavies et al.ll2004l) . 

Although there are no other reports of resolved torus emis- 
sion at this time, upper limits on t he torus size have been 
report ed in some additional sources. Prieto & Meisen heimerl 
(2004) studied a number of AGN in the 1-5 /xm range. In 
all cases the observations show unresolved nuclear emission 
at these wavelengths, setting upper limits on the torus ra- 
dius of < 5-10 pc, depending on the target distance. Even 
more significan t are t he upper limits reported at mid IR. 
Radomsk ietail d2003l) place an upper limit R < 17 pc at 10 
jim and 1 8 iim on t he nuclear component in NGC 4 1 5 1 , while 
ISoifer etaTI (|2003|) pl ace the tighter constraint R < 5 pc at 
12.5 /xm. [Soifer et al. also find an upper limit R < 14 pc for 
the 12.5 compact nuclear emission in NGC 1275. 

4.4. How Big is the Torus? 

All current observations are consistent with a torus radial 
thickness Y = R /Rd that is no more than — 20-30, and per- 
haps even as small as — 5-10. Although larger values cannot 
be ruled out, nothing in the currently available IR data re- 
quires their existence. Similarly, molecular line observations 
do not give any evidence for large toroidal structures with 
the height-to-radius ratio H/R ^ 1 required from unification 
statisti cs (see ^6*31 below). In NGC 1068. ISchinnerer et all 
( 2000) find from CO velocity dispersions that at R ~ 70 pc the 
height of the molec ular cloud distr i bution is only H — 9-10 
pc, for H/R ~ 0.15. iGalliano et all d2003l) model H 2 and CO 
emission from the same source with a clumpy molecular disk 
with radius 140 pc and scale height 20 pc, for the same H/R 
~ 0.15. Therefore, although resembling the putative torus, 
the distribution of these clouds does not meet the unification 
scheme requirement////? — 1. Evidently, the detected molec- 
ular clouds are located in a thinner disk-like structure outside 
the torus. Recent 10 iim imaging polarimetry of NGC 1068 
by iPackham et al.l d2007) shed some light on the continuity 
between the torus and the host galaxy's nuclear environments. 

As is evident from the above discussion, determining the 
torus actual endpoint is rather difficult, if not impossible; in 
fact, insisting on an endpoint for a steep 1 / r 2 distribution is 
meaningless in practice (with the currently available observa- 
tions). The torus is embedded in the central region of the host 
galaxy, and the steep radial decline of its brightness implies 
that its emission is unlikely to be cleanly separated from the 
surroundings. The only observations holding realistic chance 
for doing that are future high resolution sub-mm measure- 
ments with ALMA. Even those would require detailed anal- 
ysis that takes into account the emission from both the torus 
and its surrounding. 

It seems safe to conclude that there is no compelling evi- 
dence at this time that torus clouds beyond Y — 20-30 need 
be considered, although such large sizes cannot be excluded. 
From eq. Q] a conservative upper bound on the torus outer ra- 
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FIG. 15. — Spectral shape of the silicate 10 and 18 fim features: F\ is the torus emission in the 5-30 fim region and F c x is the smooth underlying continuum 
obtained by a spline connecting the feature-free segments in this spectral region (see text). All models have q = 2, Y = 30, a = 45°, and JVo as marked in each 
panel. Curves correspond to different tv, as labeled. Panels on the left correspond to pole-on viewing, on the right to edge-on; note the different scales of the 
vertical axes in the two cases. 



dius is then R a < 1 2 L^ 2 pc, where L45 = Lj 1 45 erg s" 1 . These 
compact dimensions have important implications for the dy- 
namics because they place the torus inside the region where 
the black hole gravity dominates over the galactic bulge. 
If the black hole mass is M.7XIO 7 M Q then it dominates 
the gravitational motions within a radius 35 (M.v/fJ 2 ) 1 / 3 pc, 
where il, typically of order unity, is the rotation velocity 
(in kms" 1 pc" 1 ) induced b y the galactic bulge in its interior 
(Elitzur & Shlosman 2006). Since the torus is well within the 
black hole sphere of gravitational influence, its dynamic ori- 
gin is determined in all likelihood by the central engine and 
its accretion disk, not by the accretion from the galaxy (see 
also jj6.6l below). 

5. SPECTRAL INDICATORS 

As noted in the Introduction, reliable analysis of the torus 
full SED in individual sources requires data that are unavail- 
able in most cases. We can expect flux measurements at wave- 
lengths longer than ~ 10 fim to be severely contaminated by 
the torus surroundings. Even in NGC 1068, the best observed 
AGN, the validity of almost all observations at A > 10 fim 
is questionable because the torus contributes only a fraction 
of the measured flux. The situation is unlikely to improve in 



the foreseeable future. The alternative approach to individ- 
ual SED fitting is analysis of large data sets in an attempt to 
identify statistical trends that might constrain the likely phys- 
ical range of torus parameters. Here we discuss the spectral 
indicators most commonly used in such analyses. 

5.1. The Silicate 10 fim Feature 

Amorphous silicate grains have strong opacity peaks due to 
the Si-O stretching and the O-Si-O bending modes, leading 
to broad features around 10 and 18 /im in the observed dust 
radiation. The 10 fim feature, the stronger of the two, is a 
common analysis tool. In smooth-density torus models the 10 
fim feature appears in emission for face-on viewing and in ab- 
sorption in edge-on viewing. As shown in §0 clumpy models 
produce more elaborate patterns. For detailed analysis of the 
features we fit a smooth curve, F c \, to the underlying contin- 
uum of the entire spectral region by a spline connecting the 
intervals 5-7 fim, 14-14.5 fim, and 25-31.5 fim. Detailed ra- 
diative transfer calculations verify that this interpolation pro- 
cedure properly reproduces the emission that would be gen- 
erated by dust stripped of its silicate features (ISirocky et al.1 
2008). Figure Q3] shows the silicate feature profiles produced 
in some representative models. As noted above, a 10 fim 
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FIG. 16. — Indicators of the 10 silicate feature (see eq.|5J: variations of the feature strength (left) and equivalent width (right) with the optical depth Ty of 
individual clouds. Model parameters are q = 2, Y = 30. Other parameters as marked. The overall optical depth along radial equatorial rays extends all the way to 
A/qtv = 4,500 in these models, yet the lOfim absorption feature is never deep. 



emission feature emerges in edge-on viewing at Ty > 100. 
This peculiarity arises because individual clouds become op- 
tically thick across the entire feature. The radiation emerging 
at this spectral range is then dominated by emission from the 
bright faces of clouds on the torus far side escaping through 
clear lines of sight. The effect becomes more pronounced as 
TVo decreases. 

The 10 /xm feature peaks at 10.0 /xm in the absorption co- 
efficients from Ossenkopf et al.l (1 19921) . and radiative transfer 
effects introduce occasional small shifts (no larger than 0.5 
/xm, mostly toward shorter wavelengths) around this value in 
the emerging spectra. To quantify the feature's strength and 
width we introduce two indicators: 



>io : 



Tn- 



EW W -- 



dX 



(3) 



7/jm 



The feature strength Sio is evaluated at the extremum near 
10.0 /xm of the continuum-subtracted spectrum. Positive val- 
ues of Sio indicate an emission feature, negative values an ab- 
sorption feature. Delineating the feature from noise in the data 
requires a certain minimum for the equivalent width EWw, 
depending on the detection system. Our sign definitions are 
matched for both indicators so that absorption produces a neg- 
ative EW\q, the opposite of the standard. 

Figure [16] displays the variations of Sio and EW\o with the 
single cloud optical depth Ty for pole-on and edge-on viewing 
and various model parameters that bracket the likely range in 
AGN tori. If the feature width were the same in all models, 
Sio and EWw would be equivalent to each other 5 , but because 
of variations in the feature shape (see figure |T3T >. £Wio con- 



If the feature is parametrized as aF c \e 



where 5 A is wavelength 



shift from the peak, then S[q= ln(l + a) while EW\q = y/naA. 



tains independent information. Figure [16] shows that pole-on 
viewing produces an emission feature only for a limited set 
of parameters. Clouds heated indirectly do not produce the 
emission feature when Ty > 20 (sec. 3.2, part I). The radia- 
tion from a directly illuminated cloud displays the feature in 
emission only toward directions with a view of a sufficiently 
large fraction of the cloud's illuminated face (see fig. 12, part 
I). An observer along the pole of a toroidal distribution will 
detect an emission feature only from direct viewing of clouds 
located within j3 < 45° from the equator. Such clouds will be 
obscured by foreground clouds in most cases, except when the 
torus width is small (er = 15°) or the overall optical depth of 
the clumpy medium is small (small Mq and Ty). Therefore, at 
i = 0°, only a small region of parameter space produces mod- 
els with a weak emission feature while most other parameters 
produce either a featureless SED or a weak absorption feature. 
It is also evident from figure[l6]that edge-on viewing is insen- 
sitive to the angular thickness a. Irrespective of optical depth, 
the absorption feature produced by a clumpy torus is never 
deep. Figure[T7]shows the variation of the two indicators with 
viewing angle for one value of Ty, a likely representative of 
actual torus clouds. 

Comparison with observ ations is hampered by the angular 
resolution problem. In the iMason et al.l d2006l) observations 
of NGC1068, the feature strength in the central 0."4, presum- 
ably dominated by the AGN torus, is Sio = -0.4. Scanning 
along the ionization cones in 0."4 steps shows large varia- 
tions in Sio and a strong asymmetry in its spatial distribution. 
Measurements with larger apertures contain significant con- 
tribution from the ionization cones, and Spitzer observations 
may be further contaminated by still larger dusty structures. 
Nevertheless, when these observations produce clear differ- 
ences between type 1 and type 2 sources, it seems reason- 
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parameters as marked. 



able to attribute such global trends to differences in viewing 
angles and to compare our model results with the observed 
trends while considering the actual numerical values only as 
guidance. The most detailed data come from the recen t com- 
pilation of Spitzer mid-IR spectra by lHao et al.l (|2007|) . Al- 
though a loosely defined sample, it is the largest gathered thus 
far, including 24 type 1 quasars, 45 Seyfert 1 and 47 Seyfert 
2 galaxies. The QSOs display almost exclusively an emis- 
sion feature with 0.45 > S\o > 0.05, but the Seyfert 1 galax- 
ies are clustered around zero feature strength, occupying the 
range 0.35 > S\o > -0.25. Almost all Seyfert 2 galaxies dis- 
play the 10 fim feature in absorption, with the distribution 
sh owing a str ong peak at -0.1 > Sw > -0.4. In addition to 
the lHao et alj results, an intriguing recent development comes 
from the Spitzer observations of seven high-luminosity type 2 
QSOs bv lSturm et all d2006l) . Although the individual spec- 
tra appear featureless, the sample average spect rum shows 
the 10 fjm feature in emission. More recently, iPolletta et al.l 
(2008) did find the feature in absorption i n a larger sam- 
ple of mid-IR selected obscured QSOs, while Weedm an et al.l 
(2006) found the 10 fjm feature either in absorption or absent 
in a sample of X-ray and mid-IR selected obscured AGN. 

A striking characteristic of all AGN spectra is the ab- 
sence of any deep 10 /um absorption features. Given the 
large optical depths implied by the X-ray data, smooth dust 
models predict very deep absorption features. Shallow ab- 
sorption features are a hallmark of clumpy dust distributions 
irrespective of geometry (part I), and the mild absorption 
strengths evident in our model results reflect this general prop- 
erty. In cont rast, ULIRGs d isplay features that reach ex- 
treme depths (lHao et al.l 12007b . This different behavior can 



be attributed to deep embeddi ng in a dust distribu tion that is 
smooth, rather than clumpy dLevenson et al.l 120071; see also 
Spoon etaLl 1200 7, Siroc kv et al.l 120081) . In principle, cold 
foreground screens intercepting the intrinsic IR emission of 
ULIRGs could also account for the deep silicate absorption 
in these sources. However, such an explanation would require 
two distinct dust components: a very optically thick dust blan- 
keting the primary radiation source and reprocessing its in- 
trinsic radiation to emerge at the enormous IR luminosities 
that identify ULIRGs, and an additional foreground screen 
that absorbs the reprocessed IR radiation to produce the deep 
silicate absorption. To remain cold, the foreground screen 
cannot provide the main reprocessing of the huge intrinsic 
luminosity, yet it must always be aligned along the line of 
sight with the primary dust blanket. Furthermore, the aligned 
screens have to be selectively associated with ULIRGs identi- 
fied with LINER- and H Il-like features because, unlike AGN, 
these sources never show shallow absorption (see fig. 1 1 in 
Siroc kvet al.ll200 8?). Such screens present a contrived solu- 
tion for the 10/im absorption in ULIRGs. In contrast, a single 
entity of smooth-density embedding dust that is both geomet- 
rically and optically thick accounts naturally for the total IR 
characteristics of deeply absorbed ULIRGs. 

Our calculations show that clumpy tori with Mo = 2 never 
produce an absorption feature and thus are ruled out for 
Seyfert galaxies, though perhaps not for quasars (see $6.4l be- 
low). The properties of the 10 /im feature found in Seyfert 
galaxies are reproduced by our models for Mo ~ 5-15, a 
~15°-45° and 7y ~30-100. When 7y increases above ~ 
100, these models produce at equatorial viewing a weak 10 
fim emission feature with a small equivalent width, offering 
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a potential explanation for the ISturm et alJ finding in QS02: 
the small equivalent width would make it hard to discern the 
feature in individual sources, bringing it out of the noise only 
in composite spectra. Therefore, if this finding is verified it 
could indicate that the optical depths of torus clouds perhaps 
are larger in QSOs than in Seyfert galaxies. However, this 
is not a unique interpretation. Another possible explanation 
is that the cloud number J\fo decreases as the luminosity in- 
creases. This point is discussed further in §16.4l below. 

5.1.1. Apparent Optical Depth 

The overall optical depth at visual along a radial ray in the 
torus equatorial plane is NqTv- With the standard dust prop- 
erties employed here, the magnitude of the optical depth at 
10/im is Tio = 0.07A/"oTy. Another quantity frequently em- 
ployed in data analysis of absorption features is the apparent 
optical depth at maximum absorption, obtained from / = e~ r " pp 
where / is the residual intensity. Therefore, from eq. [3] 
7"app,io = Sio- When the absorption is by a cold foreground 
screen that does not emit itself at these wavelengths, T app ,io 
is the actual 10/im optical depth of the screen. But when the 
absorption arises from a temperature gradient in the emitting 
dust, T app ,io can differ substantially from the actual optical 
depth, and the dependence of the two quantities on the dust 
column may bear little resemblance to each other. This is es- 
pecially true of the torus emission. As is evident from figures 
[16] and [17] the relation between tio, the actual optical depth, 
and T app ,io is multi-valued. Furthermore, although rio exceeds 
300 in these figures, T apPi io is never larger than unity. The ap- 
parent optical depth T apPj io is a poor indicator of the actual 
optical depth. 

5.2. Color Analysis 

Color-color plots, showing correlations between two col- 
ors, are a useful way to separate objects with similar 
types of spectra and r eveal underlying physical similarities. 
lAlonso-Herrero et al.l (120031) present data for nuclear fluxes 
from visual to 16 /mi for an expanded set of the CfA sam- 
ple of Seyfert galaxies. Removing all known sources of bias 



in the original CfA selection, they have constructed what is 
arguably the most complete sample of AGN currently avail- 
able. Torus observations at wavelengths up to 10 /jm are 
likely to be les s contaminated by emi s sion fr om the surround- 
ings. From the Alonso-H errero et all d2003l) data we find that 
fluxes at 1.6 /im, 3.5 /im and 10 /im provide a useful set 
of colors for comparison with our model results. Compared 
with other combinations, the models separate better with this 
choice of colors because the spectral slopes change the most 
around the selected wavelengths. Figure [18] shows colors for 
sets of torus models with a = 45°, Y = 30, q = 2 and vari- 
ous combinations of ry and No- The AGN flux is added to 
the torus emission whenever the probability for direct view 
of the nucleus exceeds 50%. In each case the colors depend 
on the viewing angle, resulting in a track of model results. 
Colors corresponding to type 1 viewing populate the upper 
right end of the track, type 2 viewing the lower left. Model 
parameters that explain the observations of the 10/im fea- 
ture also give good qualitativ e agreement with the data from 
lAlonso-Herrero et al.l (120031) . which fall inside the two re- 
gions delineated with dashed lines in the figure. While type 
2 models are spread out along the track, type 1 are grouped 
together more closely at the upper end since their spectra are 
dominated by the AGN continuum and thus are similar despite 
the broad range of parameters. 

6. ADDITIONAL IMPLICATIONS OF CLUMPINESS 

Comparison with IR observations shows that the likely 
range for optical depths of individual torus clouds is Ty ~ 
30-100 and there are A/o ~ 5-15 clouds, on average, along 
radial equatorial rays. Assuming standard dust-to-gas ratio, 
the column density of a single cloud is iV^ ~ 10 22 — 10 23 cm" 2 



and the torus equatorial column density is N x 



(eq) 



(1) 



io 2J -io 24 cm . Taking account of the torus dumpiness has 
immediate implications for a number of other issues not di- 
rectly related to its IR emission. 

6.1. The Torus Mass 
As shown in §2.3 of part I, the total mass in torus clouds 
can be written as M t0lus = muN^ jNc(r, (3) dV; note that M torus 
does not involve the volume filling factor. With the cloud dis- 
tribution from eq. [2] and taking for simplicity a sharp-edge 
angular distribution, so that the integration is analytic, the 
torus mass is M toras = 47t;«h sin a N^^R^Y I q (Y), where I q = 
1, T/(21nF) and |F for q = 2, 1 and 0, respectively. Taking 
Rd from eq.Q] the mass ratio of the torus and the central black 
hole is 

^P = 2xlO- 4 -^sina<l. 23 T/ 9 (4) 

where L^dd is the Eddington luminosity and N^ s 23 is the 

equatorial column density in 10 23 cm" 2 . Since the radial 
thickness Y is likely < 20-30 (©, the torus mass is always 
negligible in comparison with M, when q = 2. If the radial 
cloud distribution is flatter, eq.[4]may constrain the torus prop- 
erties to keep its mass below that of the black-hole. 

6.2. Total Number of Clouds 

As shown in part I, the total number of clouds, n tot , is the 
only torus property whose estimate involves the cloud size 
R c . Equivalently, R c can be replaced by the volume filling 
factor <fi, since inserting eq. [2] into eq. 3 of part I yields R c = 
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4>RJN() at the torus inner edge. If <fi is constant throughout the 
torus then n tot ~ A/"o/</> 2 for the 1/r 2 distribution, independent 
of the torus radial thickness Y. For example, if the volume 
filling factor is 10%, in order to encounter Ao = 5-10 clouds 
along each radial equatorial ray, the torus must contain n tot ~ 
10 4 -10 5 clouds, 

6.3. AGN Unification 

The classification of AGN into types 1 and 2 is based on 
the extent to which the nuclear region is visible. In its stan- 
dard formulation, the unification approach posits the viewing 
angle as the sole factor in determining the AGN type, and 
this is indeed the case for any smooth-density torus whose 
column density declines with angle (3 away from the equato- 
rial plane. The AGN is obscured from directions that have 
e -T V (/3) ^ i and v i s ibi e f rom tnose w i m e -^(P) ^ \ Be- 
cause of the steep variation of e~ T with r, the transition be- 
tween these two regions is sharp, occurring around the direc- 
tion where tv(/3) = 1. Denote this angle a then, so long as 
r v (0) > 1 and Tm{\-k) < 1, all AGN viewed at < i < \-k- a 
appear as type 1 sources, those at jir-a <i <hir as type 2. 
If f2 denotes the fraction of type 2 sources in the total popu- 
lation then f% = sin a for all smooth-density tori, irrespective 
of their specific angular profiles. This relation has been em- 
ployed in all studies of source sta tistics performed to date. 
From statistics of Seyfert galaxies ISchmitt et alj d2001l) find 
that /2 ~ 70%, hence thei r estimate a c± 4 5°. The issue is cur- 
rently unsettled because lHao et alJ (2005) have recently found 
that fi is only about 50% in Seyfert galaxies, or a ~ 30°. 

Within the clumpy torus paradigm, the difference between 
types 1 and 2 is not truly an issue of orientation but of prob- 
ability for direct view of the AGN. Since that probability 
is always finite, type 1 sources can be detected from what 
are typically considered type 2 orientations, even through the 
torus equatorial plane: if A/o = 5, for example, the prob- 
ability for that is e~ 5 = 1/148 on average. This might of- 
fer an explanation for the f ew Seyfert galaxies reported by 
lAlonso-Herrero et al.1 (120031) to show type 1 optical line spec- 
tra together with 0.4-16 /im SED that resemble type 2. Con- 
versely, if a cloud happens to obscure the AGN from an ob- 
server, that object would be classified as type 2 irrespective 
of the viewing angle. In cases of such single cloud obscu- 
ration, on occasion the cloud may move out of the line-of- 
sight, creating a clear path to the nucleus and a transition to 
type 1 spectrum. The time scale for such an event is deter- 
mined by the cloud size and velocity. Neither quantity can be 
found from the SED since optical depth is the only property 
of a single cloud that can be determined from SED analysis. 
However, at a distance r pc (in pc) from a black-hole with mass 
10 7 M. 7 (inM ), the local Keplerian speed is 208(M. 7 /r pc ) 1 / 2 
kms" 1 and resistance to tidal sheer implies that the size of 
a cloud with column density 10 23 /Vh,23 cm" 2 is restricted to 
< 10 16 Ar H) 23^ c /M. 7 (e.g. jElitzur & Snlosmanl2006h . The ra- 
tio of this cloud size and local Keplerian speed produces a 
time scale of \1 Ny^.it,^ I years, an order-of-magnitude 
estimate for a cloud crossing time across the line of sight. Al- 
though the likelihood of catching such crossing by chance is 
small, transitions between type 1 and 2 line spectra have been 
observed in a few so urces (seelAretxaga et al.ll 19991 and ref- 
erences therein), and Goodrichl (1 1 9891 [19 95) has argued that a 
couple of these cases are consistent with the change in redden- 
ing expected from cloud motion across the line of sight. It is 
worthwhile conducting monitoring observations in an attempt 



to detect additional such transitions. The most promising can- 
didates would be obscured systems with relatively small X- 
ray obscuring columns, which may minimize the number of 
clouds along the line of sight, small torus sizes, i.e., lower 
luminosities, and large black-hole masses. 
Accounting for the torus dumpiness, the fraction of type 2 

sources is fa = 1 - J^ 2 e'^ 7 ^ cos (3d (3 (eq. 9, paper I). The 
sharp-edge clumpy torus has /i = (1 -e _A/ °)sin<7, practically 
indistinguishable from a smooth-density torus when A/o ex- 
ceeds ~ 3>-A. However, the situation changes fundamentally 
for soft-edge distributions because at every viewing angle, 
the probability of obscuration increases with the number of 
clouds. As is evident from figure [T9l the Gaussian distribution 
produces a strong dependence on A/"o and significant differ- 
ences from the sharp-edge case. Since the shar p-edge angular 
distribution is ruled out by observations ( j33.U . the fraction of 
obscured sources depends not only on the torus angular width 
but also on the average number of clouds along radial rays. 
While the fraction f% = 70% requires a = 45° in the sharp-edge 
case, in a Gaussian clumpy torus it implies a = 33° when A/o = 
5 and a = 27° when A/"o = 10; in terms of the torus height and 
radius, H/R (= tana) is reduced from ~ 1 to ~ 0.7. It is note- 
worthy that the behavior of the lOfim feature in the a = 30° 
models comes closest to matching the observed averages of 
both type 1 and type 2 AGN, as is evident from figure [TBI 

6.4. A Receding Torus? 

The fraction fi of obscured AGN decreases when the bolo- 
metric luminosity increases. This has been verified in a large 
number of observations that estimate the luminosity depen- 
dence of either fi or f\ (the fraction of unobscured sources), 
or difference s between the luminosity functions of type 1 and 
2 AG N (see lHao et al.ll200l ISimpsonll2005t iMaiolino etall 
120071) . As is evident from figure[l9] the observed decrease of 
fi when L increases can be produced by either a decrease of a 
at constant A/o or a decrease of A/o at constant a. Both options 
are equally plausible because the torus inner radius increases 
as L 1 / 2 (eq. [JJ). The decreasing-cr option would arise if the 
torus height is independent of luminosity or increases more 
slowly than L 1 / 2 , the decreasing-A/o option would arise if the 
torus outer radius is independent of luminosity or increases 
more slowly than L 1 / 2 . 

The observed trend of fi to decrease with L may arise from 
either a or A/o or both. Source statistics cannot distinguish 
between the various possibilities, the only way to decide be- 
tween them is to find L-dependence in other observable quan- 
tities. The 10 /im silicate feature offers such an indicator 
( jj5.ll ). Among type 1 AGN, quasars consistently produce 
an emission feature but Seyfert galaxies are featureless on 
average, displaying either weak emission or absorption scat- 
tered around zero feature strength. In type 2 AGN the feature 
switches from clear absorption in Seyfert galaxies to appar- 
ent emission in QS02. That is, in both type 1 and type 2 
AGN the feature moves toward emission with the increase 
from Seyfert to quasar luminosities. As is evident from fig- 
ure [TBJ the decreasing-A/o option naturally produces such a 
universal trend: the feature appears in emission for both pole- 
on and edge-on viewing when A/J) decreases to ~ 2 at a fixed 
a. In contrast, the decreasing-cr option produces the observed 
trend toward stronger emission feature only in type 1 AGN, 
not in type 2; varying a has virtually no effect on the 10/iiri 
feature in type 2 viewing. Explaining the switch toward ap- 
parent emission feature in QS02 would require that in this 
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angular distribution. 



case higher luminosities not only reduce a but are also ac- 
companied by an increase in the optical depth of individual 
clouds. 

As is evident from this discussion, current observations, if 
accepted at face value and assuming that the torus contribu- 
tion dominates the 10 fim spectral range on average, can be 
explained if an increasing luminosity causes a decrease in the 
number of clouds A/o. Whether or not this is also accompa- 
nied by a decrease in the torus angular width cannot be ruled 
in or out. Obscuration statistics and the lO^tm feature do not 
yet provide decisive information to uniquely constrain the be- 
havior of the torus parameters with increasing luminosity. 

The decreasing-er scena rio is known as t he receding torus 
mo del, first suggeste d by Lawrenc3 (1199 ll) . It is intriguing 
that lArshakianl (120051) and Simpson] (120051) derived indepen- 
dently an almost identical relation tan a oc L . However, 
both studies, as well as every other analysis of obscuration 
statistics thus far, were based on sharp-edge angular obscura- 
tion. Removing this assumption affects profoundly the foun- 
dation of the receding torus model because the dependence 
on the number of clouds necessitates analysis with two free 
parameters, therefore a cannot be determined without A/o. 

6.5. X-rays and the AGN Torus 

Dusty material absorbs continuum radiation both in the 
UV/optical and X-rays, therefore the dusty torus also provides 
X-ray obscuration. But dust-free gas attenuates just the X-ray 
continuum, so clouds inside the dust sublimation radius will 
provide additional obscuration only in this band. 

Observations give overwhelming evidence for the orienta- 
tion dependent X-ray absorption expected from AGN unifi- 
cation. In general, the 2-10 keV X-ray continuum is heav- 
ily obscured in ty pe 2 sources and relatively unobscured in 
type 1 AGN (see [Maiolino & Risalitil 2007, and references 
therein). The strong orientation-dependence of the absorp- 
tion cannot be attributed to the host galaxy because the AGN 
axis, as traced by the jet position angle, is randomly ori- 
ented with respect to the galactic disk in Seyfert galaxies 
dKinnev et a l. 2000) and the nuclear dust disk in radio galax- 
ies dSchmitt et all2002l) . Yet in spite of the overall correspon- 



dence between the optical and X-ray obscuration, there is a 
significant number of AGN for which the expected charac- 
teristics are different in the two bands. Although substan- 
tial X-ray absorption is common among type 2 AGN, there 
are also unabsorbed X-ray sources that present only narrow 
emission lines in their optical spectra. Such cases can be ex- 
plained with the obser vatio nal selection effect su ggested by 
ISevergnini et all (l2003h and lSilverman et alJ (l2005h : in these 
sources, the optical light of the host galaxy outshines the AGN 
continuum and broad lin es. This suggestio n is supported by 
the su bsequent studies of Pag e et alTd2006l) and lGarcet et all 
d2007l) . The opposite case, obscuration only in X-rays, exists 
too — there are type 1 , broad line AGN with significant X-ray 
absor ption dPerola et al.l2004tlEckart et alJ2006llGarcet et alJ 
120071) . Extreme cases include quasars whose optical spectrum 
shows little or no dust extinction while their X-ray continuum 
is heavily affected by Com pton thick absorption (Brai to et alJ 
120041: iGallagher et al.f2006l) . This cannot be attributed to ob- 
servational selection effects. 

Obscuration that affects the X-rays but not the optical arises 
naturally from absorption by dust-free clouds. Conclusive ev- 
idence for such absorption comes from the short time scales 
for transit of X-ray absorbing clouds across the line of sight, 
which establish the exist ence of obscuring clouds inside the 
dust sublimation radius (Risaliti et al. 2002). Extreme cases 
involve 4 hour variability dElvis et al]|2004l) and variations in 
absorbing column of more than 10 24 cm" 2 within two days, in- 
dicating Compton thick X -ray absorption from a single cloud 
in the broad-lines region (Ris aliti et al.ll2007h . These obser- 
vations show that the torus extends inward beyond the dust 
sublimation point to some inner radius R x < R^. Clouds 
at R x < r < Ri partake in X-ray absorption but do not con- 
tribute appreciably to optical obscuration or IR emission be- 
cause they are dust-free. Since every cloud that attenuates 
the optical continuum contributes also to X-ray obscuration 
but not the other way round, the X-ray absorbing column 
always exceeds the UV/optical absorbing colum n, as ob- 
served dMaccacaro et alii 19821: iGaskell et alj|2007l) . Further, 
Maiol ino et alJ (l2001h find that the X-ray absorbing column 
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exceeds the reddening column in each member of an AGN 
sample by a factor ranging from ^3 up to ^100, implying 
that the bulk of the X-ray absorption comes from the clouds 
in t he dust-free inner port ion of the torus. This could explain 
the lGuainazzi et all J2005) finding that at least 50% of Seyfert 
2 galaxies are Compton thick. 

In steep radial distributions such as 1 / r 2 , which seems to 
adequately describe the torus dusty portion, most clouds are 
located close to the inner radius. If this radial profile contin- 
ued inward into the dust-free zone, that region would domi- 
nate the X-ray obscuration — as observed. Similar to the op- 
tical regime, the observed fraction of X-ra y absorbed AGN 
varies in versely with intrinsic luminosity (lUeda et al.l 120031: 
lHasingedl2004t lAkvlas et al.1 12006). This fraction is usually 
derived from the statistics of sources that have at least one 
X-ray obscuring cloud along the line of sight to the AGN, 
therefore it follows the behavior plotted in figure [l9]but with 
A/q corresponding to the total number of (dusty and dust-free) 
clouds. As the previous section shows, either the radial thick- 
ness a or the cloud number A/p co uld be responsible for a 
decreasing f? . iMaiolino et al.l d2007l) find that the fz fractions 
follow similar trends with L in the X-ray and optical regimes, 
indicating that whichever intrinsic parameter is responsible 
for these trends it might behave similarly in the dusty and 
dust- free portions of the torus. 

6.6. What is the Torus? 

In the ubiquitous sketch by lUrrv & Padovanil (1 1995b . the 
AGN central region, comprised of the black hole, its accre- 
tion disk and the broad-line emitting clouds, is surrounded 
by a large doughnut-like structure — the torus. This hydro- 
static object is a separate entity, presumably populated by 
molecular clouds accreted from the galaxy. Gravity controls 
the orbital motions of the clouds, but the origin of vertical 
motions capable of sustaining the "doughnut" as a hydro- 
static structure whose height is comparable to its radius was 
recognized as a problem since the first theoretical study by 
Krol ik & Begelmanl (119881) . 

Two different types of observations now show that the 
torus may be a smooth continuation of the broad lines region 
(BLR), not a separate e ntity. IR reverberation observations by 
Sugan uma et al.l d2006l) show that the dust innermost radius 
scales with luminosity as L l l 2 and is uncorrected with the 
black hole mass, demonstrating that the torus inner boundary 
is controlled by dust sublimation (eq. [TJ, not by dynamical 
processes. Moreover, in each AGN for which both data ex- 
ist, the IR time lag is the upper bound on all time lags mea- 
sured in the broad lines, a relation verified over a range of 
10 6 in luminosity. This finding shows that the BLR extends 
all the way to the inner bound ary of the dusty torus, validat- 
ing the iNetzer & Laorl (1 19931) proposal that the BLR size is 
bound ed by dust sublimati on. The other evidence is the find- 
ing by Ris aliti et al.l (120021) that the X-ray absorbing columns 
in Seyfert 2 galaxies display time variations caused by cloud 
transit across the line of sight. Most variations come from 
clouds that are dust free because of their proximity (< 0.1 
pc) to the AGN, but some involve dusty clouds at a few pc. 
Other than the different time scales for variability, there is no 
discernible difference between the dust-free and dusty X-ray 
absorbing clouds, nor are there any gaps in the distribution. 

These observations suggest that the X-ray absorption, broad 
line emission and dust obscuration and reprocessing are pro- 
duced by a single, continuous distribution of clouds. The dif- 
ferent radiative signatures merely reflect the change in cloud 



composition across the dust sublimation radius The inner 
clouds are dust free. Their gas is directly exposed to the AGN 
ionizing continuum, therefore it is atomic and ionized, pro- 
ducing the broad emission lines. The outer clouds are dusty, 
therefore their gas is shielded from the ionizing radiation, and 
the atomic line emission is quenched. Instead, these clouds 
are molecular and dusty, obscuring the optical/UV emission 
from the inner regions and emitting IR. Thus the BLR occu- 
pies r < Rd while the torus is simply the r > /?d region. Both 
regions absorb X-rays, but because most of the clouds along 
each radial ray reside in its BLR segment, that is where the 
bulk of the X-ray obscuration is produced. Since the X-ray 
obscuration region (XOR) coincides mostly with the BLR, it 
seems appropriate to name this region instead BLR/XOR. By 
the same token, since the unification torus is just the outer por- 
tion of the cloud distribution and not an independent structure, 
it is appropriate to rename it the TOR for Toroidal Obscura- 
tion Region. The close proximity of BLR and TOR clouds 
should result in cases of partial obscuration, possibly leading 
to observational constraints on cloud sizes. 

The merger of the ionized and the dusty clouds into a sin- 
gle population offers a solution to the torus vertical structure 
problem. Mounting evidence for cloud outflow (see, e.g., 
lElvisl2 004) indicates that instead of a hydrostatic "doughnut", 
the TOR is just one region in t he clumpy wind coming of f the 
black-hole accretion disk (see Elit zur & Shlosmanll2006l and 
references therein). The accretion disk appears to be fed by a 
midplane influx of cold, clumpy material from the main body 
of the galaxy. Approaching the center, conditions for devel- 
oping hydromagnetically- or radiatively-driven winds above 
this equatorial inflow become more favorable. The disk- 
wind rotating geometry provides a natural channel for angu- 
lar momentum outflow from the disk an d is found on many 
spatial scales, from protostars to AGN (Blandfo rd & Payna 
1982t lErrmiering et al.l fT99l iFerreiral 120071) . The compo- 
sition along each streamline reflects the origin of the out- 
flow material at the disk surface. The disk outer regions are 
dusty and mol ecular, as observ ed in water masers in some 
edge-on cases dGreen hill 2005). At smaller radii the dust is 
destroyed and the disk composition switches to atomic and 
ionized, producing a double-peak signature in some emis- 
sion line profiles dEracleo us 2004). The outflow from the 
atomic and ionized inner region feeds the BLR and produces 
many atomic lin e signatures, inc luding evidence for the disk 
wind geometry (lHall et al.ll2003l) . Clouds uplifted from the 
disk dusty and molecular outer region feed the TOR and 
may have been detected in water maser obser vations of Circi- 
nus (lGreenhilletal.1 120031) and NGC 3079 (iKondratko et all 
2005). Indeed. lElitzur & Shlosm an (2006) derive the cloud 
properties from constraints deduced from clumpy models for 
the IR emission and find that they provide the right conditions 
for H2O maser action. In both the inner and outer outflow re- 
gions, as the clouds rise and move away from the disk they 
expand and lose their column density, limiting the vertical 
scope of X-ray absorption, broad line emission and dust ob- 
scuration and emission. The result is a toroidal geometry for 
both the BLR/XOR and the TOR. Because of the strong pho- 
toionization heating of BLR clouds they may rise to relatively 
lower heights than the TOR dusty clouds. Detailed compar- 
isons of X-ray and optical obscuration in individual sources 
and in large samples should help to constrain the parameters 
A/o, fx and t v separately for the TOR and the BLR/XOR. Such 
comparisons must consider the large scatter of obscuration in 
individual sources around the sample mean (see paper I, §4.2). 
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In the outflow scenario, the TOR disappears when the bolo- 
metric luminosity decreases below ~ 10 42 erg s" 1 because the 
accretion onto the central b lack hole can no longer sustain the 
requir ed cloud outflow rate (El itzur & Shlosmanl2 006; Elitzur 
2007). With further luminosity decrease, suppression of cloud 
outflow spreads radial ly inward a nd the BLR, too, disappears. 
The recent review bv lHol d2008l) presents extensive observa- 
tional evidence for the disappearance of the torus and the BLR 
in low luminosity AGN. 

The Circinus Seyfert 2 core provides the best glimpse of the 
AGN dusty/molecular component. Wa ter masers trace bot h 
a Keplerian disk and a disk outflow dGreenhill et al.ll2003b . 
Dust emission at 8-13/im shows a disk embedd ed in a slightly 
coole r and larger, geometrically thick torus dTristram et alj 
2007). The dusty disk coincides with the maser disk in both 
orientation and size. The outflow masers trace only parts 
of the torus. The lack of full coverage can be attributed to 
the selectivity of maser operation — strong emission requires 
both pump action to invert the maser molecules in individual 
clouds and coincidence along the l ine of sight in both position 
and velocity of two maser clouds dKartie et al.lll999l) . Proper 
motion measurements and comparisons of the disk and out- 
flow masers offer a most promising means to probe the struc- 
ture and motion of TOR clouds. 

7. SUMMARY AND DISCUSSION 

We have developed a formalism for handling radiative 
transfer in clumpy media and applied it to the IR emission 
from the AGN dusty torus. In the calculations we execute 
only the first two steps of the full iteration procedure out- 
lined in §3.2, paper I, and the moderate total number of clouds 
considered here validates this procedure. When that number 
increases, the probability for unhindered view of the AGN 
decreases, the role of indirectly heated clouds becomes more 
prominent and eventually requires higher order iterations . Our 
current calculations employ some additional simplifying ap- 
proximations: The grain mixture is handled in the composite- 
grain approximation, all dust is in clouds without an inter- 
cloud medium and all clouds are identical. We have already 
begun work on removing these assumptions and will report 
the results in future publications. 

In contrast with the smooth-density case, the clumpy prob- 
lem is not well defined because clouds can have arbitrary 
shapes, and any given set of parameters can have many 
individual realizations. Our formalism invokes a statisti- 
cal approach for calculating an average behavior, and it is 
encouraging that other approache s produce similar results. 
iDullemond & van Bemmelf (120051) conduct "quasi-clumpy" 
calculations in which the torus is modeled as a set of axisym- 
metric rings, and compare the results with the smooth-density 
case. In agreement with our conclusions they find that only 
smooth-density models can produce very deep absorption fea- 
ture while clumpy dust produces stronger near- IR, broader 
SEP a nd much more isotropic IR emission. Hfani g et ail 
(2006) employ 3D Monte carlo calculations that bypass some 
of our approximations. They also treat different cloud realiza- 
tions for the same global parameters, allowing them to show 
the intrinsic scatter in SED due to the stochastic nature of the 
problem. Their results are in agreement with ours, validat- 
ing our approach and the approximations we employ. Since 
the du st properties in their calculations are from lDraine & Lea 
(1984), the 10 /im feature reaches somewhat larger strengths 
than in our calculations, which employ the Ossenko pf et alj 
d 19921) "cool" dust (but are similar to our original results 



in iNenkova et al. I [20021 which also employed iDraine & Led 
dust). In spite of these differences, Honig eTalJ (l2006h too 
find that the silicate absorption feature is never as deep as ex- 
pected for a uniform dust distribution, and obtain qualitatively 
similar behavior of the silicate emission feature and overall 
SED shape. 

The models presented here show that clumpy torus models 
are consistent with current AGN observations if they contain 
No ~5-15 dusty clouds along radial equatorial rays, each with 
an optical depth ry ~30-100. The cloud angular distribution 
should decline smoothly toward the axis, for example, a Gaus- 
sian profile centered on the equatorial plane. Power-law radial 
distributions r~ l - r~ 2 produce adequate results. Dust grains 
with optical properties of the standard Galactic mixture pro- 
vide satisfactory explanation to the IR observations. The be- 
havior of the 10/im silicate feature, in particular the lack of 
any deep absorption features, is reproduced naturally without 
the need to invoke any special dust properties. Several sugges- 
tions that the abundance or composition of AGN dust might 
differ from its Galactic counte rpart can be discard ed because 
of subsequent developments. iRisaliti et alj d 1999b note that, 
assuming standard dust abundance, the large column densities 
discovered in X-ray absorption imply torus masses in excess 
of the dynamical mass, posing a problem for the system stabil- 
ity. However, their mass estimates were based on the uniform 
mass distribution and large torus sizes derived from smooth- 
density models. The compact sizes and steep density distri - 
butions of clumpy m odels eliminate the problem (see 36.11 ). 
Maiol ino et alJ (1200 ll) suggested that the widely different UV 
and X-ray extinctions they found in individual sources could 
imply low dust abundance, but the subsequent discovery of 
rapid variations shows that X-ray obscuration by dust-free 
clouds is the more likely explanation (see 36.51 ). They also 
invoked the lack of prominent 10/im absorption features as an 
indication that AGN dust is different from Galactic, but this is 
a natural consequence of clumpy dust distributions (see 35.11 ). 
Intrinsic extinction curves deduce d from spectral analysis of 
type 1 sources (see ICzerny||2007l for a recent review and a 
comprehensive discussion of uncertainties) generally indicate 
a depletion of small grains, as could be expected: the obscu- 
ration in type 1 sources is dominated by the dusty clouds clos- 
est to the center and these clouds contain predominantly large 
grains, which survive at the smallest distances from the AGN 
(see 32.11 ). There is no compelling evidence for significant 
differences between the properties of AGN and Galactic dust. 
Other dust compositions are not ruled out, but nothing in the 
current data requires major departures from the dust grains we 
use. 

The close proximity of dust temperatures as different as > 
800 K and ^200-300 K found in interferometry around 12/xm 
cannot be explained by smooth-density models even when 
they account fo r the individual tempe ratures of grains with 
different sizes (Schartmann et al. 2005). Clumpiness resolves 
this puzzling observation because the dust on the dark side 
of an optically thick cloud is much cooler than on the bright 
side. Thanks to the mixture of different dust temperatures at 
the same radial distance, clumpy models naturally explain the 
torus compact size. In spite of the high anisotropy of its ob- 
scuration, the torus emission is observed to be nearly isotropic 
at A > 12 /im. Clumpy models resolve this puzzle too, since 
the emission from a torus with radial thickness Y = 10 varies 
little with viewing angle. The variation is especially small if 
the radial distribution is 1 /r 2 or steeper, and such steep radial 
profiles maintain a nearly isotropic emission even at larger 



AGN Dusty Tori: II. Observational Implications 



19 



torus sizes. 

In addition to IR observations, dumpiness significantly im- 
pacts the analysis of other data, in particular obscuration 
statistics. The fraction fi of obscured sources is controlled 
not only by the torus angular thickness a, as in all analyses to 
date, but also by the cloud number Afo- With Ao = 5, a 70% 
fraction of type 2 AGN implies a ~ 30° instead of the stan- 
dard 45°. Observations indicate that increasing the bolometric 
luminosity from the Seyfert to the quasar regime induces (1) 
a decrease of fi and (2) a switch to emission feature at 10 
/im for both type 1 and some type 2 AGN. Both trends can 
be explained with a change in a single torus parameter — A/q 
decreases from ~ 5 in Seyfert galaxies to ~ 2 in QSO (see 
figures IT61 and [T9l). Decreasing a, the scenario known as the 
receding torus model, explains the first trend but has no ef- 
fect on the second. The emergence of the 10 /J,m in emission 
would require in this case the additional increase of individual 
clouds optical depth to Ty > 100 in QSO. 

The decreasing-A/o scenario provides the simplest explana- 
tion for the trends observed when L is increasing, but that 
does not guarantee its validity. This demonstrates the diffi- 
culties in deducing the model parameters from observations 
that cannot yet resolve the torus basic ingredient — the indi- 
vidual dusty clouds. The problem is compounded by the lack 
of angular resolution that hinders clean separation of the torus 
component from the flux measured at most IR wavelengths 
and by the degeneracies of the radiative transfer solutions that 
prevent decisive, one-to-one associations between model pa- 
rameters and observable quantities. The only practical way 
around these difficulties is to match trends identified in the 
data with similar general properties of the models. 

Our main conclusions can be summarized as follows: 

• The torus angular distribution has to be soft edged 

• Clumpy models can produce nearly isotropic IR emis- 
sion together with extremely anisotropic obscuration 

• Clumpy models can explain all current observations 
with compact torus sizes; SED fitting is a poor con- 
straint on the size 



• Standard interstellar dust describes adequately AGN 
observations; there does not seem to be a need for any 
major modifications of grain properties 

• Clumpy sources never produce a very deep silicate fea- 
ture; apparent optical depth, obtained from / = e" T "pp 
where / is the residual intensity at maximum absorp- 
tion, is a poor indicator of the actual optical depth 

• The probability for direct line-of-sight to the AGN at 
large viewing angles is small, but not zero 

• The statistics of obscured sources depend on both the 
torus angular thickness and the number of clouds along 
radial rays 

• The torus and the BLR are the dusty (outer) and dust- 
free (inner) regions in a continuous cloud distribution; 
a more appropriate designation for the torus is Toroidal 
Obscuration Region (TOR) 

• X-ray obscuration comes from both TOR and, predom- 
inantly, BLR clouds 

As long as IR observations are incapable of resolving indi- 
vidual torus clouds, one must rely on the combined evidence 
for clumpy structure instead of on a "smoking gun". Individ- 
ual TOR clouds seem to have been resolved in observations 
of outflow water masers in Circinus and NGC 3079. Proper 
motion measurements and comparison of these masers with 
their disk counterparts provide the most promising method 
for probing the TOR structure and kinematics. The Circinus 
AGN, whose dust emission has been resolved in VLTI obser- 
vations, is an especially attractive target for studying the dusty 
and molecular content of TOR clouds. 
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APPENDIX 
TECHNICALITIES 

The relevant coordinates in describing both the cloud distribution and the source function are the cloud's radial distance r, 
angle (3 from the equatorial plane and the angle a between its radius vector and the AGN-observer axis (see eq. [2] and part I 
figure 2 and eq. 8). The torus emission requires an integration along a path inclined by the viewing angle ; from the torus axis at 
some displacement from the center (eq. 5, part I). To handle the geometry we introduce a cartesian coordinate system centered 
on the AGN with z toward the observer and x—y in the plane of the sky, with the torus axis in the y-z plane at angle i from the 
z axis. The integration path is specified by its fixed values of x and y, so that the angular displacement is (8 x ,8 y ) = (x/D,y/D) 

and the angular impact parameter in brightness profiles is 8 = ^/x 2 + y 2 /D. The integration variable is z. At any point r = (x,y,z) 
along the path, the cloud coordinates are found from 



tan/3 = ■ 



ysini+zcosi 



, cosa=-. (Al) 

yx 2 + Cycos«-zsin/) 2 r 

The path integration in eq. 5, part I, produces the intensity generated by the cloud distribution. Since our source function 
calculations involve only the first two steps of the full iteration procedure described in §3.2, part I, we must introduce a correction 
to take proper account of flux conservation. With pagn the fraction of the AGN luminosity that gets through the torus (eq. 8, part 
I), we calculate I^(x,y;i), the brightness map of clumpy torus emission in the direction i, from 



L(l-p ACN ) 



where H x (x,y;i) = P eS c.\(r)S Ct x(F)N c (r)dz 



(A2) 



4tt J dcosi J dX J H\dxdy 

Here S Ci \(r) and Nc(r) are, respectively, the source function and column density of clouds at position r along the integration path, 
and P e sc,\(r) is the probability for a photon of frequency A to escape from that point through the rest of the path. The torus flux 
at distance D and viewing angle i is calculated from F^(i) = (1 /D 2 ) J I^(x,y; i)dxdy. With these expressions, the spectral shape is 
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determined from the first two iteration steps while ensuring that the torus emission properly obeys flux conservation (eq. 17, part 
I). 

The quantity H\ is intrinsically a function of scaled variables, H\ = H\(x/ Rj,y /R&i) , because the brightnes s at position (x,y) 
depends only on the distribution of dust in temperature and optical depth along the path (llvezic & Elitzur|[l997l) . Therefore, from 
eq. IA2l the brightness has the form I\(0 X , 6 y ; i) = (L/47r/?d 2 )/(#<:/#d,#y/#dX where / is a dimensionless function of the scaled 
angular displacements. Since the brightness scale L/'XitRa 1 is determined by the dust sublimation temperature r su b (eq.[T} , the 
only effect of the luminosity is to set the overall angular scale 9^, effecting a self-similar stretch of the brightness map. Similarly, 
the flux, F^, is a product of the bolometric flux Fagn and a luminosity-independent spectral shape. 
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